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Abstract 

QJ ' Multi scales method is used to analyze a nonlinear differential-difference 

^^ I equation. In the order e'^ the NLS eq. is found to determine the space-time 

evolution of the leading amplitude. In the next order this has to satisfy a 
complex niKdV eq. (the next in the NLS hierarchy) in order to eliminate 
secular terms. The zero dispersion point case is also analyzed and the 
relevant equation is a modified NLS eq. with a third order derivative term 

K^ ' included. 

^^ I Many one-dimensional systems of biological interest are very complicated 

^—N ' structures, formed from complexes of atoms - we shall call them "molecules" 

- connected by hydrogen bounds. It is usually assumed that only one of the 

2^ . intra-molecular excitations plays an active role in the storage and transport 

of energy in these systems. In the case of a-helix structure in protein this 
corresponds to the amide I vibration (C=0 stretching). We shall call this intra- 
molecular excitation the vibronic field. Localized excitations of solitonic type 
can exist in these systems, due to a nonlinear interaction between the vibronic 
O ■ field and the acoustic phonon field, describing the molecule oscillations along 

the chain. The simplest model starts from a Frolich Hamiltonian and with an 
ansatz - coherent state approximation - for the state vector describing this type 

p\ ' of localized excitation ( Q - [^ and references therein) . 

j^ ■ After eliminating phonon variables a nonlinear differential (time) - difference 

(space) equation for the vibronic coordinate is obtained; 

L{{yn}) = G({y„}) (1) 

where L is the linear part and G the nonlinear one. For the specific example 
we have in mind, originating from Takeno's model P], L and G are given by 

HiVn}) = -^ + ^IVn - Y^ JmnVm (2) 

ra^n 



a^ 



G{{yn]) = ^AyrAvl+i + yl-i) - Bvl (3) 



1 

2' 

In the linear part the last term is a long range interaction between vibrons, and 
we shall assume that Jmn decreases exponentially (Kac-Baker model) 



J^n = Jlm-n\ = ^L V""""'"""' ' ' = """• ^^^ 
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The first term in r.h.s. of G results from the nonlinear interaction between 
vibrons and phonons, while the second one from a quartic anharmonicity in the 
vibron Hamiltonian. 

The linear equation admit plane wave solutions e*^, 9 = kan — uit, where 
u){k) is given by the dispersion relation 

D{Lok) = LO (k) — \ loq — 2 2_. Jp cos kap 
V P=i 

2 f .2 .2 ^-"^ ( sinh7 



o;^ - Un - ^Ir^— ' , - 1 (5) 

V ° ^^ 2r Vcosh7-cosA;a )) ^ ' 

describing an optical vibrational branch with uj'^(k) a monotonously increasing 
function of k from tJ^(O) = uJq —oj\^ to ^^^(vr) = ujq +uj\^j^. We shall assume 
that a no-resonance condition takes place 

Dy = D{vLO, vk) / 0, u(^N\ v^l. (6) 

It is well known that the effect of a weak nonlinearity occurs at large space- 
time scales, determining a redistribution of energy on higher harmonics, and a 
modulation of amplitude. In order to investigate these effects we shall use the 
multi-scales method (reductive perturbation method) Q. The method starts 
by introducing slow space-time variables 

X = ean, tj = eH (7) 

and expanding y„ in an asymptotic perturbative series, which due to the form 
(3) of the nonlinearity G is given by 

odd oo 

yn = J2 e"' E e'^P,-(^; *1' ^2, ...) + c.c (8) 

u=l p=v 

Recently several papers have used this method to discuss the propagation of 
quasi- monochromatic waves in weakly nonlinear media l5|-[0], or of long surface 
waves in shallow waters Q . Very interesting are the conclusions concerning the 
role played by the NLS hierarchy , or the KdV hierarchy Q in eliminating the 
secular terms which would destroy the asymptotic character of the perturbative 
series. Of special interest for the present paper is the reference P], which will 
be followed as close as possible. 

In calculating the time derivative we have to take into account that t appears 
in as well as in the slow time variables ii, i2v Also in writing the expressions 
for yn±i,ym we have to expand the corresponding amplitudes around the point 
n. Taken these precautions the calculations are straightforward (although quite 
tedious in the higher orders): the asymptotic expansion is introduced in (1), 
(2), (3) and the coefficient of each power of e and each harmonic e^^^ is equated 
with zero. 

In the first order in e we re-obtain the dispersion relation (5). In the order 
e^ the amplitude Yi^i has to satisfy the equation 



and consequently Yi^i will depend only on the variable ^ = x — vt, where 
Vg = -^ = oji is the group velocity. 

In the next order e^, from the terms proportional with e*^ we get 



L+Y2,i = —^ - i^2(n,i) (10) 




K2{Yi,i) = ico2 -^^ + q\Yi,i\'Yi^i (11) 



Here u;„ = ^^ and q = ^(2 + cos2ka - 3f ). As the r.h.s. of (10) is in 
the null space of L-|-, l2,i will blow up linearly in ti unless the r.h.s. is strictly 
equal with zero, i.e. Yi^i has to evolve in t2 according to the cubic nonlinear 
Schrodinger equation (c = y~) 



Oto 



(^ + 2c|y,,|Xi) (12) 



In this case 1^2,1 will depend also on the characteristic coordinate ^ only. From 
terms proportional with the third harmonic e one obtains 

^3^3,3 + {A cos ?,ka - B)Yl^ = (13) 

and due to the no-resonance condition (6) it is an algebraic equation giving 
13^3 in terms of ^1,1. The same thing happens for all the higher harmonics 
and the corresponding amplitudes Yp^y can be explicitly written in terms of 
Ypi and their derivatives. Therefore we shall concentrate our attention to the 
amplitudes ^,1, related to the first harmonic e . 

The solution of the NLS eq. (12) depends on the sign of LO2 and q. As u)i 
vanishes at /c = and k = it, there is a point kc £ (0, vr) for which uj2{kc) = 0. If 
k < kc{k > kc) we have uJ2 > ^{^2 < 0). The sign of q depends on the constants 
A and B. For A > and i? < -j it is always positive, while for S > ^ it is 
negative. Depending on the sign of uj2 and q the NLS eq. (12) can have bright 
or dark soliton solutions. 

In the order e^ from the terms proportional with e we get 

9Y2,i j^,,^ . ayi,i a^yi.i . wiu;2^ d\Yi^ ^dYi^i 

Ot2 Ots Ot,-^ UJ Ot, O^ 

(14) 
Here 

K'2{Yi,i) = iuj [-^ + MYlxYl^ + 2|yi,i|2y2,i) j (15) 

is the Frechet derivative of K2, and q\ = -^- The l.h.s. of (14) is the linearized 
NLS eq. It is well known that the commuting symmetries cJj of the NLS eq. 
are solutions of this equation. As they are important for our further discussion 
we remained the expression of the first ones (by ^ we shall denote a solution 
of the NLS eq.) 

do = 

cT2 = i(^ + 2c|*|^*) cj3 = -(^ + 6c|^|^^) (16) 



i^. 


d^ 
"■=8? 


t;:+-i*p*) 


-3--(5^3+6c|^l%^ 



The eq. (14) is a forced linear equation for ^2,1- It is necessary to identify 
secular terms in the r.h.s. of (14) and then to fix the is dependence of li^i 
in such a way to eliminate their effect. These secular terms have to be found 
between the members of the null space of linearized NLS eq., i.e. between 
the commuting symmetries ctj. Indeed if such a symmetry a would exists it 
will generate a t2(T contribution to ^,1, and the asymptotic character of the 
expansion (8) would be destroyed in a time ^2 = 0{e~^). Two such symmetries 
((To,(T3) are easily seen in the r.h.s. of (14), if it is written in the form 

- ^ + ^3 (^ + 6c\Y,,,\'Y,,^^ + N{Y,,,) (17) 

where 

iOli02^ 5|Yi,i|2 ,2^^1,1 



N{Y,,,) = -2c^^Y,,^^^l^ + {q, - 6ca;3)|^i.i 



In order to avoid this secular behaviour we require that the ts dependence of 
Yi^i is given by the following complex modified KdV equation 



dt: 



f + -3(^ + 6c|y,,|X,)=o 



which is the next equation in the NLS hierarchy. The influence of the rest 
A^(li^i) on l2,i can be further treated using a Green function formalism |^]. 
Let us consider a single soliton solution flC 



M/ = 2-i^ (19) 

y c cosh z 

(/>(e, t2) = 2SiC + 4u;2(S? - P?)t2 + <Po 

z{C,t2) = 2Pi{^-C0 + AuJ2Slt2) 

where Si , Pi are the real and imaginary part of the complex eigenvalue Ci = 
Si + iPi in the inverse scattering transform method, and </>o,^o ^-re the initial 
phase and the initial position of the soliton. Applying the above procedure, 
in order to eliminate the possible secularities, the soliton parameters must be 
ts-dependent. This dependence can be found introducing (19) in (18). The 
complex eigenvalue will remain unchanged, while for i?i'0;?o the following linear 
equations are found ||^] 

^ = -8^3^1(5? -3Pf) 

^ = -8uJsPi{:iSf - P?) (20) 



A similar analysis was given by Kodama |11], and the same results are obtained 
using the direct perturbation method of Keener and McLaughlin |^. More 
complex situations and details will be published elsewhere. 



Let us consider now the situation when u;2 = 0, i.e. the propagation of a 
wave with the wave vector kc- As loi has a maximum at this point it represents 
the wave propagating with the highest group velocity. A similar situation is 
encountered in the case of pulse propagation in nonlinear optical fibers where 
this point is known as the "zero dispersion point" (ZDP) |jl2|-|l^]. The power 
required to generate an optical soliton is minimal in this point, and its evolution 
in space and time is governed by a modified NLS eq., with a third order deriva- 
tive included. We shall show that a similar situation appears in the present 
case. 

In applying the multiple scale method we shall use the same asymptotic 
expansion (8) for the vibronic variable. Then the nonlinearity contribution 
begins with terms of order e^. To have contributions of the same order from 
third order derivatives we have to change the scaling of the S, variable, namely 

^ = el (an -Wit). (21) 

We have to take into account also a dependence of the phase 6 of the propagating 
wave on the slow variable ^. Defining the local wave number k as the derivative 
of the phase 9 with respect to (an) we find that k is slightly different from kc, 
and the simplest choice is 

k = kc{l + e^). (22) 

Expanding all the quantities depending on k around the point kc in the order 
e^ the following equation is found for the leading amplitude Yi^i — > ^, 

i^T + 3^xx + i^xxx + (3|^P^ = (23) 

where X = kc^, T = Qt2, ^ = oj^k^ and Q = ^. It has the same form as 
the equation describing the propagation of nonlinear pulses in optical fibers in 
the ZDP region ||l2|| -[p^]. In our case it makes the transition between the two 
regions, where bright and dark solitons exist. It seems that it is not completely 
integrable, but some analytical and numerical results suggest that some long- 
living localized excitations exist [0]. Further investigations are necessary. 
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